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ABSTRACT 


Contact stress distribution of spiral bevel gears 
using non-linear finite element . static analysis is presented. 
Procedures have been developed to solve the non- linear 
equations that identify the gear and pinion surface 
coordinates based on the kinematics of the cutting process; 
and, orientate the pinion and the gear in space to mesh with 
each other. Contact is simulated by connecting GAP elements 
along the intersection of a line from each pinion point 
(parallel to the normal at the contact point) with the gear 
surface. A three - dimensional model with four gear teeth 
and three pinion teeth is used to determine the contact 
stresses at two different contact positions in a spiral bevel 
gearset. A summary of the elliptical contact stress 
distribution is given. This information will be helpful to 
helicopter and aircraft transmission designers who need to 
minimize weight of the transmission and maximize 
reliability. 
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CHAPTER I 
INTRODUCTION 

Spiral bevel gears are used to transmit power between 
intersecting shafts. One critical application is a 
helicopter transmission in which power is transmitted from 
the horizontal engine to the vertical rotor shaft. 

These gears operate at relatively high rotational 
speeds and transmit substantial power, (i. e. 15,000 rpm and 
900 hp) Accurate prediction of tooth bending and contact 
stresses are important to design stronger, lighter gears with 
longer service lives. 

Much effort has focused on predicting stresses in 
gears with the finite element method. Most of this work has 
involved parallel axis gears with two dimensional models. 
Parallel axis components have closed form solutions that 
determine surface coordinates. Spiral bevel gears, however, 
do not have closed form solutions available to identify 
surface coordinates. 

Only a few researchers have investigated finite 
element analysis of spiral bevel gears. [1-3] The research 
reported here will develop a finite element model to study 
meshing of a spiral bevel gearset. Pinion and gear tooth 
surfaces will be developed based on the kinematics of gear 
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manufacture. The individual teeth are then rotated in space 
to create a multi -tooth model. 

The initial model has one pinion tooth and one gear 
tooth. Subsequent models consist of four gear teeth and 
three pinion teeth. The tooth pair contact zones are modeled 
with GAP elements. The model development procedures and 
finite element results are presented. 
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CHAPTER II 


GEAR SURFACE GEOMETRY 

This chapter describes the gear manufacturing 
process, the kinematics of cutting, tooth surface coordinate 
solution procedure, concave and convex surface rotations of 
the gear and pinion, and different orientations required for 
the spiral bevel gears to mesh with each other. Also 
described is the transformation process which transforms a 
point on the cutting blade to a point on the workpiece. 
Application of tooth surface coordinate solution technique is 
also described. 

2.1 GEAR MANUFACTURE 

The Gleason Works, Rochester, NY, provides machinery 
for the manufacture of spiral bevel gears. The machine shown 
in figure (l) is used to cut these gears [4] . These machines 
are preferred because they can be used for both milling and 
grinding operations. Grinding is especially important for 
producing hardened high quality aircraft gears. This machine 
consists of three main parts: the machine frame, the cradle 
and the sliding base. The gear cutter is mounted to the 
cradle of the cutting machine. The machine cradle with the 
cutter may be imagined as a crown gear that meshes with the 
gear being cut as shown in figure (2) . The gear which is 
being cut and the cradle with the mounted head cutter rotates 
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slowly about its own axis. The combined process generates the 
gear tooth surface. 

The cradle rotates far enough so that only one space 
between the tooth is cut out. It then reverses rapidly while 
the workpiece is withdrawn from the cutter and indexed ahead 
by the translation of the sliding base (on which the work is 
mounted) with respect to the cutter, in preparation for the 
next tooth. This sequence of operations is repeated till the 
last tooth is cut [4] . 

The head cutter, which holds the cutting blades or 
the grinding wheels, is shown in figure (2) . The cutter 
carries blades of straight lined- prof ile which generate a 
cone during rotation. It rotates about its own axis C-C at a 
speed for efficient metal removal, independent of the cradle 
or workpiece rotations. The head cutter is connected to the 
cradle through an eccentric that allows adjustment of the 
axial distance between the cutter center and the machine 
center. The adjustment of the angular position between the 
two axes provides the desired spiral angle. 

The pinion is typically cut on one side at a time, concave 
and convex tooth sides, whereas the gear is cut both sides 
simultaneously (duplex method) . The cradle and workpiece are 
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connected through a system of gears and shafts, which 
controls the ratio of rotational motion between the two. This 
ratio is called the ratio of roll. For cutting, the ratio is 
constant; but for grinding it is a variable. 


2.2 TOOTH SURFACE COORDINATE SOLUTION PROCEDURE 

From the kinematics of the manufacturing process an 
analytical method was developed by Litvin [4,5] to generate 
the surface coordinates of spiral bevel gears. The 
manufacturing data, the surface geometry and the equation of 
meshing are used to develop the equations. 

The first equation is the EQUATION OF MESHING given 
in reference [4-6] . This equation is based on the kinematics 
of manufacture and the machine tool settings. The equation of 
meshing requires that the relative velocity between a point 
on the cutting surface and the same point on the pinion being 
cut must be perpendicular to the cutting surface normal. 

This is given by: 

n • V = 0 

where 'n 1 is the normal vector to the cutter and the 
workpiece surfaces at the specified location of interest, and 
•V 1 is the relative velocity between the cutter and the 
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workpiece surfaces at the specified location. This equation 
is developed in terms of the machine tool coordinates u,0,4>c 
where u and 0 (length and rotational orientation) locate a 
point on the cutting head cone and <J>c is the rotated 

orientation of the cutter as it swings on the cradle. This 
equation which is in terms of (u, 0, <J>c ) could be transformed 

into (x, y, z) coordinate system using the transformation 
process described in the next section. The equation of 
meshing for straight sided cutters with a constant ratio of 
roll and left-hand spiral tooth is given by [4-6] 

(u - rcotip cosip) cosy sinx 
+ s[(m cw - siny)cosip sin0 

- cosy sinip sin(q - 4>c ) 3 (2.2.1) 

+ E m (cosy sinip + siny cosip cost) 

- iijn siny cosip sinx = 0 

Where y is the root angle of the gear being manufactured, ip 
the blade angle, q the cradle angle and x equals (0 - q + 

<J>c) • 

Since the equation of meshing has three unknown 
parameters u, 0, <j) c , two more equations have to be developed 

to solve for points on the surfaces. The two other equations 
are obtained considering the pinion and gear geometry. 

The orientation of the gear on a coordinate system with 
the apex of the cone at the origin is as shown in figure (3) . 
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Any point on this surface is located by determining 
'r' and 'z' coordinates when projected on to the XY plane. 
These are determined using the geometry of the gear. 

The value of Z W/ found by the transformation process 
of a point from the head cutter to the workpiece surface, 
should match the value of Zbar found by the projection of the ’ 
tooth in the XY plane. To satisfy this condition the 
following equation is used. [1] 

Z - Zbar = 0 2.2.2) 

The radial location of a point from the work axis of 
rotation must be equal to the value of r as shown in figure 
3. This forms the third equation [1] 

R bar ” V ( Xw + Yw ) .......( 2.2.3) 


These three above equations could be written as 

fi(u, 0, <|> c ) = 0 
f2 (u, 0 , <J) C ) = 0 

f 3 (u , 0 , 4 > c ) = 0 (2.2.4) 

2 . 3 COORDINATE TRANSFORMATIONS 

. The purpose of the coordinate transformations are to 
transform a point on the cutting head to a point on the 
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workpiece. These points on the workpiece will then describe 
the gear surface geometry. During the cutting process, the 
cutting head acts as if in mesh with a simulated crown gear. 
Therefore, the appropriate point in the cutting head is a 
point that satisfies the equation of meshing. (i. e. the 
normal of the point on the generating surface should be 
perpendicular to the relative velocity between the cutter and 
the gear tooth surface.) 

To perform these transformations homogeneous 
coordinates are used to allow rotations and translations of 
vectors by simply multiplying the matrix transformations. To 
transform a point on the cutting head to the workpiece five 
such transformations are carried out. These are as follows 
[4-6] : 


1. Cutter coordinates 

(S c ) 

--> 

Cradle coordinates 

(S s ) 

2. Cradle coordinates 

(S s ) 

--> 

Machine fixed 


(Sa) 

3 .Machine fixed 

(Sm) 

--> 

Pinion-pitch angle 

(S p ) 

4 .Pinion-pitch angle 

(S p ) 

--> 

Pinion axis of 

rotation 

(S a ) 

5. Pinion axis of rotation 

(S a ) - 

-> Pinion fixed 

system 

(S w ) 


Before discussing these the following points should 
be noted: 

1. The head cutter coordinate system S c is rigidly connected 
to the coordinate system S s . The system S s is rigidly 
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connected to the cradle that rotates about the X m axis of the 
machine coordinate system s m . 

2. The cutting head is attached to the cradle with a specific 
orientation 's' and ’q' where 's' is the distance between s a 
and s c origins (s = 0 C 0 S ) and <j> c is the roll angle of the 
cradle. 

3 . The cradle rotates at a constant (0 (co cra cii e or oo c ) . 

4 . The workpiece rotates at a constant (0 (toworkpiece or tow) . 

5 . The ratio of to cra die/ ^workpiece is the same as the ratio of 

angular velocity of the simulated crown gear meshing with the 
gear being cut. 

6. The cutting head rotation is not relevant. It only effects 
the metal removal rates. It is helpful to think of the 
cutting head as being fixed. The cradle rotation sweeps the 
cutting head cone through the workpiece cone. The 
interference between the two shapes the tooth surface. 

7 . The orientations of the coordinate transformations are 
different for left-hand and right-hand gears. The Gear tooth 
spirals to the left while looking from the front of the gear 
(the front is being viewed from the apex) in case of left- 
hand gear whereas, for the right-hand gear it spirals to the 
right (see figure 4) . In the following discussions the left- 
hand gear transformations are described exclusively. 

Beginning with a point on the cutting blade given by 

[4-6] : 
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r c 


rcos tp - ucos ip c * 
usinipcsin 6 
usin ip c cos 0 
1 


(2.3.1) 


where 'r ' is the radius of the blade at x c =0 in figure (2), 
and 'ip ' is the blade angle. Parameters u and 0 locate a point 
in the system S c and are unknowns, whose values will be 
determined. 


To transform the head cutter into coordinate system 
S s , the following matrix is used : 



i 

o 

o 

o 


o 

cosq 

sinq 

0 


0 

-sinq 

cosq 

0 


0 ' 
-s sinq 
scosq 

1 


(2.3.2) 


where ' q 1 is the cradle angle and 1 s 1 is the distance between 
the coordinate systems S s and S c as shown in figure 5. In 
this view, x m the axis of rotation of the cradle, and x c the 
axis of rotation of the cutting head are coming out of the 
page. These are parallel to each other and fixed in space. 

's' and 'q' are polar coordinates in the S s coordinate system 
that locate the cutting head. They are different for each 
gear design. Once the cutting head is located on the cradle 
for a specific design, it is locked in place with 's' and 'q' 
being fixed and constant for the manufacture of that gear. 
[Msc] * ic is the coordinate of a point in the S s coordinate 
system. 
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To transform S s to the fixed coordinate system Sm the 
following transformation is used: 


M 


ms 


‘10 0 O' 

0 cos <|> c sin <j> c 0 
0 -sin<(> c cos<j) c 0 

0 0 0 1 


(2.3.3) 


where <j> c is the roll angle of the cradle. The coordinate 
system S m is fixed in space. The origin of Sm is at the 
center of rotation of the cradle. The origins O m and O a 
coincide as shown in the figure (5) . The rotation of Ss 
relative to S m indicates how much the cradle has rotated. 
[Mms] [Msc] * f c are the coordinates of a point in the 
coordinate system S m . 


The third coordinate transformation is from the 
coordinate system S m to S p which orients the pitch apex of 
the gear being manufactured. This transformation requires 
special machine tool settings 1^ and E m along with the 
dedendum angle which is obtained from the design data. The 
values of L m and E n settings could be obtained from table I. 
The coordinate system orientation to generate a left hand 
gear surface is shown in figure (5) . 

The transformation matrix is given by: 
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M 


pm 


cos 8 0 -sin 8 
0 1 0 
sin 8 0 cos 8 
0 0 0 


-L m sin8‘ 
+E m 

, cos 8 


m 


1 


(2.3.4) 


The next transformation rotates a point from S p to 
S a . This is well illustrated in figure 6. The common origin 
of S p and S a locates the apex of the gear under consideration 
with respect to the S m coordinate system. This requires a 
rotation about Y a by the pitch angle |i . This is given by the 
matrix shown below: 



cos |i 0 -sin [i 0 
0 1 0 0 

-sinn 0 cosjx 0 

0 0 0 1 


(2.3.5) 


The final transformation is from the coordinate 
system S a to coordinate system S w , which is fixed to the 
component being manufactured. A rotation about the Z p axis is 
required through an angle <Jw» as shown in figure 6. This 
transformation is given by: 


M 


wa 


'cos<t> w 

- sin<j> w 


0 

0 


sin <|> w 0 O' 

cos <|> w 0 0 

0 1 0 

0 0 1 


(2.3.6) 


Using these five transformations any point in the head cutter 
coordinate system can be transformed into the workpiece 
coordinate system with this equation [4-6] : 
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A w 


IMJ [Map] [Mp m ] [M ms ] [Ms C ] 


rcos tp - ucos ipc 
usin ipcsin 0 
usin ip c cos 0 
1 


(2.3.7) 


2.4 APPLICATION OP SOLUTION TECHNIQUE 

The three equations discussed earlier to describe the 
surface coordinates are non-linear equations which do not 
have a closed form solution. These are solved using a 
numerical technique. These three equations could be rewritten 
as 


F (u, 0, 4> c ) 


fi (u, 0, <i> c ) =0 

f 2 (u / 0, 4>c) = 0 

_f 3 (u / 0, 4> c ) = 0_ 


(2.4.1) 


or 


F(x) = 0 

where 


(2.4.2) 



(2.4.3) 


The equation of meshing which is written in terms of 
( u, -0, <J> C ), is numerically differentiated by adding small 

increments using a five point formula [7], and fhon 
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coordinate transformations are used to convert the change in 
u, 0, <|> c into change in x,y,z. By doing so the values of 


afi 0fi flfi 
dxi ’ 0X2 ’5x3 

are calculated. The other two equations are already in the 
x,y,z coordinates. Newtons method [7] is then used to solve 
the equations . The procedure is as follows: 

Starting with an initial guess x° for the three parameters 
Xi, X2 and X 3 / F(x°) is calculated. From this we obtain 

■fl(X V 

F(X°) = f 2 (X°) (2.4.4) 

-f 3 (X °) . 


Now, calculating the Jacobian for X°, nine specific values 
are obtained. 

af j (x °) af ] (x °) af j (x °) 

5x x 5X2 5 x 3 

5f 2 (X °) 5f 2 (X °) 8 f 2 (X °) 

8xi ax 2 5x 2 

5f 3 (X °) 5f 3 (X °) 5f 3 (X °) 

5xjl 5x 2 5x 3 

This Jacobian is formed by numerical differentiation and 
evaluating the derivatives at the initial guess. 

The set of non-linear equations can now be 
transformed into a set of linear equations. This is done as 
follows: 
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at, (x °) 

0fi(x °) 

0f 1 (X °) 

axj_ 

0X 2 

0X 3 

0f 2 (X °) 

0f 2 (X °) 

0f 2 (x °) 

dXi 

w 

(0 

ax 3 

0f 3 (x °) 

0f 3 (X °) 

af 3 (x °) 

dXx 

0X 2 

ax 3 


f i (X °) 


f 2 (X °) 


I- Y 3 _j L f 3 (X °) j 

The solution to this set of linear equations is used 
as the guess to obtain the solution in the next iteration. 
This process is repeated until the difference between the 
solutions in two consecutive iterations is within the 
required tolerance. 


x 1 = x° + y 

or; u 1 = u° + Yi 

e 1 = e ° + y 2 
<t> c 1= 4>c° + y 3 

The variables used to define the above parameters are 
given in reference [3] which illustrates the numerical 
example solved for a specific design data. 


2.5 CONCAVE AND CONVEX ORIENTATIONS 

The orientations are performed in order to obtain 
the top land on the gear and pinion tooth. The top land 
thickness is calculated at the face angle on the toe end of 
the gear tooth. In figure 7, PI is the concave side location 
of the face angle point at the toe end of the tooth. 
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P2 and P2' are the initial and final convex locations of the 
face angle point at the toe end of the tooth [3] . These 
points are described in coordinate system S w . T a and Tb are 
the desired and initial top land thickness for the gear 
tooth. By rotating the convex side points of the gear tooth 
by an angle e the desired thickness for the gear tooth is 

obtained. 

2.6 GEAR AND PINION ORIENTATIONS 


The apex of the pinion and the gear cones generated 
meet at a point as shown in figure 8. To place the gear and 
pinion in mesh with each other the following rotations are 
carried out. 


(i) The pinion is rotated about its own axis by a very small 
rotation angle in the counter clockwise direction about the 
z axis. The amount of rotation is determined by observing 
contact with the three dimensional geometric modeling program 
(TDGMP) [8] . The amount of rotation used for the example 
presented in Chapter V is 3.56 degrees. 

(ii) The second transformation is the rotation of all the 
surface points obtained from the previous transformation 
about the global Y axis by 90 degrees in the counter 
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clockwise direction. This transformation places the pinion 
over the gear as shown in figure 8. 

(iii) The third transformation is to rotate the gear tooth 
surface points by. (360/Nt) + 180 degrees about the global Z 
axis in the clockwise direction. By doing so the gear tooth 
generated, meshes with the pinion as shown. The rotation 
used in this study was 190 degrees. 
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CHAPTER III 


CONTACT SIMULATION WITH GAP ELEMENTS 

After the gear and the pinion are oriented in mesh 
with each other, the contact is simulated by connecting GAP 
elements. An assumed contact point is found on the pinion 
concave surface and then the normal vector is calculated from 
this point. The intersection of the normal vector with the 
gear convex surface is found. The intersection of parallel 
vectors from neighboring grid points are also found. These 
points are connected to form the GAP elements. The solution 
technique which performs these steps is presented in this 
chapter. 

3.1 PROCEDURE TO FIND THE CONTACT POINT 

After orienting the pinion concave surface to be in 
mesh with the gear convex surface, the two surfaces are 
viewed in TDGMP [8] to find an approximate area of contact. 
This is done at various slices across the section of the gear 
surfaces in the direction of the spiral angle as shown in the 
figure 9. Having found the area of contact, the two closest 
grid points are chosen to be the contact point. 

3.2 PROCEDURE TO DETERMINE THE NORMAL VECTOR 

The equations used to find out the normal vector is 
provided by Litvin [4]. The normal vector is first found in 
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the head- cutter coordinate system and then is transformed to 
the workpiece coordinate system using the transformations 
described in chapter II. 

Any point in the tool surface which is a cone could 
be represented in the coordinate system S c which is rigidly 
connected to the head-cutter is as follows: 

rcos ip - ucos ipc 
usin ip c sin 0 
usin tp c cos 0 

1 J (3.2.1) 

The coordinate system S s is an auxiliary coordinate 
system that is also rigidly connected to the tool as 
mentioned earlier. The following matrix transformation is 
used to represent the generating surface in coordinate system 
Ss (for a left hand gear) . 

'1 0 0 0 iK' 

0 cosq -sin q -ssinq y c 
0 sin q cos q s sin q z 
0 0 0 ! J[lJ ( 3 . 2 . 2 ) 




, Equations (3.2.1) and (3.2.2) yields 

x s = r c cot'?,, - u cos^ c 

y s = u sin^ c sin(0 - q) - s sinq 
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= u sinip c cos(0 - q) + s cosq 


(3.2.3) 


z 


The unit normal to the generating surface in coordinate 
system S s given by [6] 


n s = 


i^r where 



dfc 

5u 


(3.2.4) 


From equation (3.2.4): 


— = Oi + usintp c cos (0 -q) j - u sinip c sin ( 0 - q) k 
50 (3.2.5) 


^2. = -cos ip c i + sin ip c sin (0 - q) j + sin n> c cos( 0 - q) k 
5u (3.2.6) 


using (3.2.2 ) to ( 3.2.6): 


N. = 


l 

0 

-COS ip c 


u sin tpc cos(0 - q) 
sinipc sin(0 - q) 


-u simpc sin(0 - q) 
sinip c cos(0 - q) 


(3.2.7) 


or 


N7 = [ u sin 2 ip c cos 2 (0 - q) + u sin 2 ip c sin 2 (0 - q) ]i + 

[ u simpc cos ipc sin ( 0 - q) cos ( 0 - q) ] j + 

[ u sin ipc cos ipc cos ( 0 - q ) ] k 


(3.2.8) 


and 

pSfsF = u 2 [ sin 4 ip c + sin 2 ip c cos 2 ip 0 sin 2 (0 - q) + sin 2 v c cos 2 ip c cos 2 (0 - q)] 
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or 


N s = Vu 2 sin 2 ip c = u sin tp c 


( 3 . 2 . 9 ) 


Therefore; 


n s = & = sin ipc i + cosip c [ sin (0 - q) j + cos (0 - q) k] 
PM 


or 


n s 


sin ip c 

cos tp c sin ( 0 - q ) 
cosipcCos (0 -q) 


( 3 . 2 . 10 ) 


In order to define the unit normal with reference to 
the machine fixed coordinate system S m t the following 
procedure is used. First any surface point with respect to 
S m coordinate system is given by: 



'X s ' 

y s 

LZsJ 


( 3 . 2 . 11 ) 


where 




0 

cos<(> c 
—sin <> c 
0 


0 

sin<t» c 

cosdi 


( 3 . 2 . 12 ) 


The unit normal with respect to S m coordinate system is given 


as: 

[flm] = Pms] [ n s] 


( 3 . 2 . 13 ) 
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where [L mg ] may be determined by deleting the fourth column 
and row in matrix (3.2.12). Thus 



1 

0 

0 

sin V c 

n m = 

0 

cos<(> c 

sin <j) c 

cos \y c sin (9-q) 


0 

-sin <t> c 

COS(|> c _ 

cos\y c cos(6-q) 


(3.2.14) 


which gives 




sin ip c 
cos \p c sin t 
cos tp c cos T 


(3.2.15) 


where t = 0 - q + <j> c for the pinion which is a left hand 
gear. Note that the normal vector is a function of 0 and <j) c 
since q, the cradle angle, is a constant for a given point. 

In figure 10, [n] is seen to be the same along the cutting 
head cone for any u. 

Equation (3.2.15) gives the normal vector with 
reference to the machine fixed coordinate system. In order to 
generate it with respect to the workpiece coordinate system 
the following transformation is carried out: 


n w — [Mpm] [Map] 


sin ip c 
cos ip c sin t 

COS \p c COS T 
1 


(3.2.16) 
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3.3 PROCEDURE TO DETERMINE THE INTERSECTION OF NORMAL 
VECTOR WITH GEAR SURFACE 

After finding the normal vector with respect to the 
workpiece coordinate system it is required to find the 
intersection of the normal vector (from the pinion concave 
surface) with the gear convex surface. To do so the pinion 
concave surface (and the normal vector at the contact 
point) is oriented to be in mesh with the gear convex 
surface. During the procedure the gear convex surface is 
fixed because the equations to solve for the intersection 
point are in terms of u, 0 , <j)c in S w . The pinion is rotated 

to the gear. The different rotations carried out on the 
pinion are described below and are shown in figure 11. 

STEP 1: The pinion concave surface, with attached normal 
vector, are rotated by an angle such that the pinion and gear 
do not interfere during meshing. (Note: the algorithm 
produces the gear and pinion with random rotational 
orientation) 

STEP 2: The pinion concave surface and attached normal vector 
are rotated by 90 degrees in counter clockwise direction 
about the Y axis 
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STEP 3: The pinion concave surface, with attached normal, are 
rotated by 190 degrees counter clockwise about the z axis to 
location 3 in figure 10 (to be near the gear convex 
surface) . (Note that meshing actually occurs at location 2 . ) 
The rotation in step 3 to location 3 is for convenience in 
solving for the intersection point . This location is 
convenient because .u, 0 , 4>c identify a gear surface point in 

S w at location 3. 

STEP 4: The final rotation of the pinion concave surface with 
attached normal is a slight counter clockwise rotation about 
the Z axis. This is done to compensate for the rotation of 
the gear convex surface done to create the top land. 

The procedure used to find the intersecting point is as 
follows: 

Consider the coordinates of the contact point to be 
(X pin , Y p i n , Z pin ) as shown in figure (12) . Let a point in 

space which the normal vector passes through be given as 
(Q X / Qy/ Qz) and the unit normal vector coefficients be 

given as (n x , n y , n z ) . Then: 

( Qx > Qy 9 Qz ) = ( Xpin 9 Ypin , Zpin ) + ^ ( n x , Ily , Iljr) (3.3.1) 
where 'b' is a scalar factor. 
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Scale factor b is found by incrementing its value (see figure 
12) such that it becomes a point on the gear convex surface. 
For this to happen it has to satisfy the following 
equations given below. The first equation is given as : 


Qx = Xgear = Xpj n + b (n x ) 

Qy = Ygear = ^pin + (ily) 

Qz = Zgear = Zpin + b (ilz) (3.3.2) 

Any point in the gear convex surface could be found by using 
the relation as given below (provided the machine settings 
are maintained for this surface generation) : 


Xgear 



rcos ip - ucos tpc 

Ygear 

II 

L 

[Map] [Mpm] [Mms] [Ms C ] 

usin ip c sin 0 

•fgear 



usin tp c cos 0 

1 J 



1 


Let 

[Mwa] [Map] [Mpm] \Mms] [M sc 3 


then using (3.3.2) and ( 3.3.3) 


"Xpin + b (n x ) 

Ypin + b (ily) 

^pin + b (II2) 
1 


m 


rcos ip - ucos ipc" 
usin ip c sin 0 
usin ip c cos 0 
. . 1 


( 3.3.4) 


or 


M 


'rcos ip - ucos ipc 


"Xpin + b (n x ) 

usin ip c sin 6 


Ypin + b (lly) 

usin tpcCOS 6 


Zpin + b (n^ 

1 


1 J 


= o 


(3.3.5) 


From equation (3.3.5) three equations could be obtained such 
that 


Ai(<l>,b,u) = A 2 (<J>,b,u,0) = A 3 (<g>,b,u,6) = 0 (3.3.6) 

By using the equation of meshing given in equation (2.2.1) 
along with the three other equations given in (3.3.6), the 
unknowns u, 0, 4)c, and b can now be solved for. These values 
of u, 6 , <frc, b will now define the intersection of the normal 
vector with the gear convex surface. The same procedure is 
adopted to find the intersection of all vectors (parallel to 
the contact point normal, and from grid points on the pinion 
concave surface) with the gear convex surface. Later, these 
intersecting points are rotated in TDGMP [8] so as to lie on 
the gear convex surface in mesh. The following rotations are 
carried out: 

STEP 1: The intersection points are rotated by an angle 

required to create the top land of the gear surface 
about Z axis. 
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STEP 2: Later these points are rotated by (360/N c ) + 180 

degrees in the counter clockwise direction to lie on 
the surface of gear convex surface in mesh. 

The intersecting points which do not lie within the boundary 
of the gear convex surface are discarded. Note that while 
doing this care should be taken so as not to create excessive 
distortion of the hex elements on the gear. 

3.4 VERIFICATION OF THE INTERSECTION POINT 

To check the results of the procedure adopted earlier, 
the following method is used. The (u, 0, <j) c ) values obtained 

at the intersection of the gear convex surface is transformed 
into (Xgear/ Ygear, Z gea r) values using the coordinate 
transformation given below. Note that initial machine 
parameters for the gear convex surface generation should be 
used during this transformation. This is given by equation 
(3.3.3) which is: 


Xg^ 


rcos ip - ucos ipc 

Ygear 

— [Mwa] [Map] [Mpm] [Mns] [Msc] 

usin ip c sin 6 

Zgear 

1 J 

usin ipecos 0 
1 


The result obtained is now substituted back in 
equation (3.3.2) 
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Xgcar = Xpi n + b (n x ) 

Ygcar = Ypjn + b (lly) 

Zgear = Zpin + b faz) 

The value of ( Xgear/ "^gea.t, ^gear ) obtained/ if 
same, verifies the point at which the normal vector pierces 
the gear convex surface and is the intersection point of 
interest . , 

The above method is used to locate the intersection of 
normals from all points on the pinion surface (in the contact 
zone) with the gear surface. The gear tooth surface is 
remeshed utilizing the intersection points as shown in figure 
13. GAP elements are connected between corresponding nodal 
points on the pinion and the intersection points on the gear 
surface. 
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CHAPTER IV 


NUMERICAL EXAMPLES 

This chapter describes a numerical example to determine 
the values of the contact stresses of spiral bevel gearset 
using an example design. Three different models of these 
gears at different levels of difficulty are analyzed. The 
first model simulates a two tooth contact, one pinion and one 
gear tooth. The second model simulates two contact points 
using four gear teeth and three pinion teeth. The third model 
simulates the two contact points along with added hub 
effects. Various contact stress levels are discussed. 

4.1 MODEL DESCRIPTIONS 

All models were constructed in a TDGMP [8] from the 
output of the computer programs and were analyzed for contact 
stress distribution using a finite element program [9] . A 
non-linear static analysis solution technique with GAP 
elements is used for the analysis. The modeling was done on 
a workstation and the analysis was done via a supercomputer . 
The models where constructed using 8 noded HEX elements and 2 
noded GAP elements . 
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MODEL I: The first model contains one pinion and one gear 
tooth as shown in figure (14) . The model consists of 1280 
elements and 1980 nodes. (5940 degrees of freedom) . The 
torque load was applied as a concentrated force in the beam 
as shown. In all of the models; the pinion was fixed by 
specifying zero displacement on the four comer nodes of each 
rim face, 8 nodes total were fixed in this manner; and the 
gear was free to rotate about its axis of rotation (z axis) . 

Model I contains 15 GAP elements along the contact 
area. Initially, only one is closed. This simulates the 
contact condition in the model. The beam section has a very 
high modulus to prevent excessive bending. The various stress 
results obtained are as shown in figures (15-18) . The 
solution iterates 5 times to confirm the assumption on the 
GAP elements (i.e. open or closed) is in agreement with the 
finite element displacement solution. The structure of the 
contact area and the history of the GAP elements at each 
iteration is given in figure (19) . 

The values of the minimum element principal stresses as 
plotted by the TDGMP [8] and the average values of the 
minimum principal stresses at various nodes in the contact 
area is shown in table (I) . 

MODEL II: This model contains four gear teeth and three 
pinion teeth in mesh. The model consists of 7626 nodes and 
5292 hex elements. (22878 degrees of freedom) . The contact 
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stresses at two different contact regions are analized with 
this model. A metal huh is provided on the gear to transmit 
the forces to the two different contact areas. The model 
generated and the boundary conditions imposed are as shown in 
figure (20) . The model contains 15 GAP elements in the major 
contact region and 6 on the edge contact region. ' 

The GAP elements in the major contact region are 
located using the algorithms described in this report. At 
the edge contact, the GAP elements are located by plotting 
lines from the pinion surface and finding the intersection 
point with the gear surface with the TDGMP [8] . In both 
contact regions one GAP element is initially closed and the 
rest are open. The solution iterated 4 times to reach 
equilibrium. The status of the GAP elements after every 
iteration is shown in figure (29) . The values of the min im u m 
principal stresses and the average minimum nodal principal 
stresses are given in table (II) . The stress contours are 
also shown in figures (21-28) . 

MODEL III: This model contains four gear teeth and three 
pinion teeth, similar to Model II, with the added hub region 
as shown in figure (30) . It contains 10 101 nodes and 7596 
elements (30303 degrees of freedom) The model contains 15 
GAP elements in the first contact region and 6 in the edge 
contact. As in MODEL II one GAP element in both the contact 
regions is closed initially, and the others open. The 
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solution iterated 3 times before obtaining the final solution on 
the fourth iteration. Different stress contours are shown in 
figures 31-34. The status of each GAP element after each iteration 
is shown in figure 35. Table III gives the element minimum 
principal stress values and the average nodal minimum principal 
stress values. 

For all three models the torque is applied as a force at a 
moment arm of 2.0976 inches. Hence, the total torque applied on 
the gear is approximately 9911 inch pounds. The pinion torque is 
three times less or 3303 inch pounds. 
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CHAPTER V 


DISCUSSION OF RESULTS 


MODEL I: This model consists of one gear tooth and one 
pinion tooth with minimal hub. In this model the element 
minimum principal contact stress on the pinion is found to be 
around -195,923 psi. This value occurs at the node connected 
by GAP element 1231 as shown in figure (19) . The average 
nodal minimum principal stresses showed a higher value of 
about 380,000 psi, and with a maximum of 447,700 psi.. The 
maximum influence of average nodal minimum principal stress 
is at the node connecting GAP element 1228. Six GAP elements 
are closed in the final iteration. 

The nodal stresses are higher because they occur at the 
GAP elements which acts like a concentrated force applied 
directly to a node. The elemental stresses are based on 
nodal averages and, as a result, are lower than the nodal 
stresses. 

Model II: This model was created to distribute the load to a 
second pair of teeth in contact. The minimum element 
principal stresses at the major contact region was found to 
be -156,041 psi. The average nodal minimum principal stress 
is found to be around -300,000 psi. The maximum influence is 
at the node connected by GAP element 1232 (see figure 29) 
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which showed a value of -339,064 psi. The stresses 
decrease, compared to Model I, because of load distribution 
to the edge contact region. 

At the edge contact, the minimum principal element 
stress is found to be around -101,661 psi and the average 
nodal stresses are found to be about -220,000 psi. The 
maximum average nodal minimum principal stress value was 
concentrated at the node connected by GAP element 6200 (see 
figure 29) is around -435,002 psi. 

The GAP element status at each iteration shows that 
initially 9 GAP element at the major contact area and 1 at 
the edge contact are closed. Also at the time of 
convergence, 5 in the major contact area and l in the edge 
contact area are closed. The' 5 GAP elements closed were also 
closed in the final iteration of model I. 

Model XXI: This model also consists of four gear teeth and 
three pinion teeth in mesh. Also modeled is a bigger section 
of hub. (see figure 30) The minimum principal element stress 
obtained from the TDGMP [8] is found to be -123,967 psi. The 
mean nodal principal stresses obtained at the contact region 
is found to be 260,000 psi. The maximum influence is at node 
connecting, GAP element 1232 (see figure 35) and the value 
is -299,931 psi. At the edge contact region, the minimum 
principal element stress is found to be around -84,932 psi. 
The mean of the average nodal minimum principal stresses in 
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the contact region is found to be -200,000 psi. The maximum 
influence is at node connected by GAP element 6200 ( see 
figure 35) and it is -363,720 psi. These results are 
consistent with estimates of hertzian contact stresses for 
these gears. 

At the end of the first iteration, 8 GAP elements in 
the major contact area and 1 in the edge contact area are 
closed and after the final iteration 4 in the major contact 
area and l in the edge contact area are closed. The same GAP 
elements closed in the major contact region were also closed 
in Models I and II. 

For this model the area of the contact region is 
marked as shown in figure 36. To estimate the area the 
average area due to the influence of each GAP element is 
taken as shown in figure (37) . The area of each rectangle is 
found to be 0.0049 sq. in. Hence, the area of the entire 
contact region is 4 x 0.0049 which equals 0.0196 square 
inch. in* in. 

For all three models, large node to node stress 
variation indicates a need for mesh refinement. 
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CHAPTER VI 


CONCLUSIONS 

The authors presented a method for doing nonlinear finite 
element contact stress analysis of spiral bevel gears. The 
method incorporates the following features. 

1) A model was developed for obtaining spiral bevel gear 
surface points. This model is based on numerical analysis of 
the manufacturing kinematics. 

2) The gear and pinion are orientated for mesh with a series 
of coordinate transformations . 

3) The gear surface is remeshed to provide proper 
orientation of the GAP elements used to simulate contact. 

4) A series of models were analyzed with increased 
complexity.- 

A) Model I - One pinion and one gear tooth in mesh, 

(5940 degrees of freedom) 

B) Model II - Four gear teeth and three pinion teeth in 
mesh. (22878 degrees of freedom) 

C) Model II - Four gear teeth and three pinon teeth in 
mesh with added hub region. (30303 degrees of freedom) 

5) Results agreed well with calculated Hertzian contact 
stresses. 

6) Large node to node stress variation indicate a need for 
refinement. 
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Figure 3 - Projection of gear tooth into XY plane 
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4. Right and left hand spiral bevel gears. 
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Figure 8. - Gear and pinion orientations required for meshing. 
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Figure 11. - Rotations required to solve for intersection of pinion 
normals with gear surface. Location 1: Small rotation about z axis 
(prevents interference during meshing). Location 2: -90 degree 
rotation about Y axis. Location 3: -190 degree rotation about Z 
axis to place pinion in mesh with gear. Also; small' rotation of 
gear convex surface about Z axis to create top land. 
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from pinion intersecting gear 
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Figure 13 . - Gear geometry distortion after connecting gap 
elements. 


50 


Figure 14 . Model I ? pinion and gear tooth in mesh with 

boundary conditions. 
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Figure 15. Model I; minimum principal element stresses 
in pinion tooth. 
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Figure 16. - Model I; Minimum principal element stresses 

in pinion. 
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Figure 17. - Model I; average nodal minimum principal stresses in 
pinion tooth. 


54 


Psi 

A -18707. 


B -48295. 


-77883. 



D -107471. 

£ -137060. 

F -166648. 
G -196236. 
H -225824. 
X -255413. 
j -285001. 
K -314589. 


Figure 18 - - Model I ; average nodal minimum principal stresses in 
pinion at contact region. 
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FIGURE 19.- Model I: GAP element status at each iteration. 
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Figure 20. - Model II. Mesh generation with boundary conditions. 
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Figure 21. - Model II. Minimum principal element stresses in 
pinion having major contact. 
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Figure 22. Model IX; minimum principal element stresses at maior 
contact region of pinion tooth. 
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23. Model XI. Average nodal minimum principal stresses in 
pinion having major contact 
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FicfU3r e > 24 . Model XX. Average nodal minimum principal stresses in 
pinion having edge contact . 
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Figure 25. - Model II. Minimum principal element stresses in pinion 
having edge contact. 
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Figure 26.- Model II. Minimum principal element stresses 
at edge contact region in pinion. 
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Figure 27. - Average nodal minimum prinicipal stresses in pinion 
having edge contact. 
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Figure 28 . - Model II. Average nodal minimum principal stresses in 
pinion at edge contact region. 
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FIGURE 29. Model II: GAP element status at each iteration. 
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Figure 30. - Model III. Finite element model and boundary 
conditions 
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Figure 31. - Model III. Element minimum principal stresses in 
pinion at major contact region. 
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Figure 32. - Model III. Nodal minimum principal stresses in 
pinion at major contact region. 
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Figure 33. - Model III. Element principal stresses in pinion 
having edge contact 
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Figure 34. - Model III. Nodal principal stresses in pinion 
having edge contact 
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FIGURE 35.- Model III: GAP element status at each iteration. 
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DISTANCE BETWEEN GRIDS 56-46 
DISTANCE BETWEEN GRIDS 46-47 
APPROX. AREA OF EACH DASHED BOX 
TOTAL AREA OF CONTACT 


= 0.14302 in 
= 0.035118 in 

= 0.14032 * 0.035118 = 0.00492 sq.in 
= 4 • 0.00492 = 0.01968 sq.in 


• GAP ELEMENTS CLOSED AFTER FINAL ITERATION 
AREA OF INFLUENCE OF GAP ELEMENTS 


Figure 37 - Contact area calculation 
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TABLE I 


MA.TOR CONTACT AREA STRESSES 


AVERAGE NODAL MINIMUM PRINCIPAL STRESSES (Psi) 


- 99,278 

- 447,732 

- 292,902 

- 128,307 

- 32,075 

- 73,283 

- 113,703 

- 330,360 

- 292,452 

- 87, 131 

- 67,939 

- 94,100 

- 99,650 

- 58,195 

- 209,021 


MINIMUM PRINCIPAL ELEMENT STRESS VALUES (Psi) 


- 76,077 

■■EH 

- 195,923 

- 120,920 



inn 

BHKH 

mm 


- 58,335 

1HRE9 

- 98,768 

- 85,134 

- 59,703 


Each box represents a grid, point where the gap element is connected 
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TABLE II 


MAJOR CONTACT AREA STRESSES 

AVERAGE NODAL MINIMUM PRINCIPAL STRESSES (Psi) 




^BSHHB 

- 80,869 

iww 

- 47,977 


BMW 

- 170,588 

- 22,499 



BHHB 

- 35,836 



MINIMUM PRINCIPAL ELEMENT STRESS VALUES (Psi) 



- 104,243 

- 156,041 

- 97,885 

- 21,234 

- 38,595 

HU 

- 148,742 

- 120,439 

- 29,393 

bbbb 

- 51,781 

- 77,700 

- 62,839 

- 20,692 


ROOR CONTACT AREA STRESSES 

AVERAGE NODAL MINIMUM PRINCIPAL STRESSES (PSI) 


BBBB 

- 47,887 

BUM 

- 123,496 

- 435,002 

- 481,34 


MINIMUM PRINCIPAL ELEMENT STRESSES (PSI) 


- 25,826 

- 67,699 

Bum 


- 101, 661 

- 42,817 


Each box represents a grid, point where the GAP element is connected 
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TABLE III 


MAJOR CONTACT AREA STRESSES 


AVERAGE NODAL MINIMUM PRINCIPAL STRESSES (Psi) 


- 39,532 


HHI 

- 51,752 

- 17,504 

■HU 

HHi 

Hi 

iRKHHi 

- 16,417 

- 32,710 

HEH 



- 9,714 


MINIMUM PRINCIPAL ELEMENT STRESS VALUES (Psi) 


HlH 

Hi 

- 123,967 

- 75,732 

- 13,044 

- 29,618 

- 60,895 

- 120,019 

- 90,893 

- 17,253 

- 28,952 

- 41,204 

- 62,573 

- 47,993 

- 140,90 


EDGE CONTACT AREA STRESSES 

AVERAGE NODAL MINIMUM PRINCIPAL STRESSES (PSI) . 


HH 


- 185,999 


- 363,720 

- 39,929 


MINIMUM PRINCIPAL ELEMENT STRESSES (PSI) 




HH 





Each box represents a grid point where the gap element is connected 
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